Interface free energies in p-spin glass models 
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The replica method has been used to calculate the interface free energy associated with the 
change from periodic to anti-periodic boundary conditions in finite-dimensional p-spin glass models 
in the phase which at mean-field level has one-step replica symmetry breaking (IRSB). In any finite 
dimension the interface free energy is exponentially small for a large system. This result implies 
that in finite dimensions, the IRSB state does not exist, as it is destroyed by thermal excitation of 
arbitrarily large droplets. The implications of this for the theory of structural glasses are discussed. 
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The use of p-spin glass models as models for under- 
standing structural glasses was pioneered in the work of 
Kirkpatrick, Thirumalai and Wolynes They found 
that for the infinite-range version of these models, for 
which mean- field theory is exact, that there were two im- 
portant temperatures. At a temperature there exists 
a dynamical transition below which an ergodicity tran- 
sition occurs. The dynamical equations which arise are 
very similar to those of the mode-coupling theory of liq- 
uids 0- At a lower temperature Tc there is a first-order 
(discontinuous) transition to a phase with a IRSB order 
parameter. As T Tc the entropy goes to zero, so this 
transition is identified with the Kauzmann glass transi- 
tion Tk the temperature at which the extrapolated 
entropy of the supercooled liquid apparently falls to that 
of the crystalline phase. 

It has been realized for many years that outside the 
mean-field limit there will be no genuine dynamical tran- 
sition at Td- This is because the metastable states which 
trap the system dynamically are unstable because of acti- 
vation processes over the finite free energy barriers which 
exist outside the mean-field limit. However, the appar- 
ent divergence of the viscosity in fragile glasses at some 
non-zero temperature, such as that in the Vogel-Fulcher 
formula Ijl 
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has encouraged the belief that there might be a real ther- 
modynamic glass transition near Typ- The observation 
that the ratio Tk/Tvf lies beteen 0.9-1.1 for many glass 
formers whose Tk ranges from 50 K to 1000 K [5| is 
further evidence for this possibility. However, we have 
recently argued that there is no genuine glass transition 
and that Tj^ and Typ are just 'crossover tempera- 
tures'. The Ising spin glass in a field seems to provide a 
satisfactory analogy for supercooled liquids near Typ. 

We shall here strengthen these arguments by showing 
that the IRSB state simply does not exist in the finite 
dimensional p-spin glass model! The reasons are similar 
to those why there is also no real dynamical transition 
Td in finite dimensional systems. There is no transition 



at Tc to a state with one-step replica symmetry breaking 
because that state is destabilized by thermal fluctuations 
within it of droplets of arbitrarily large size which only 
cost a finite amount of free energy to produce. When 
p is even, the Hamiltonian is left unchanged by fiipping 
the signs of all the spins, and the destablizing droplets 
are simple to identify: They are droplets of the time- 
reversed state. Our basic task then is to calculate the 
interface free energy for droplets. 

There have recently been other attempts to calculate 
interface free energies in p-spin models [3, 1^ ■ Their re- 
sults are very different to those reported here, essentially 
because they are studying a different interface free en- 
ergy - that appropriate for non-equilibrium effects below 
Td- The interface free energy studied here is the one com- 
monly studied in Ising spin glasses which involves chang- 
ing the boundary conditions along one direction (the z 
direction) from periodic to anti-periodic 0. That is, we 

study SF = \J AFp (here and in the following, the 
overbar means averaging over bond configurations) where 
AFp AP — Fp — Fap, and Fp and Fap are the free ener- 
gies with periodic and anti-periodic boundary conditions 
respectively. Anti-periodic boundary conditions can be 
realized by reversing the sign of the bonds crossing a 
plane whose normal is parallel to the given direction. It 
will be shown below that AFp ^p = 0. If the system is of 
length L in the z direction, we find that the interface free 
energy SF decreases to zero as ~ exp(— L/2f), provided 
L ^ 1^. £, is the correlation length in the IRSB state. 
This state will therefore be destabilized by the thermal 
excitation of large droplets, since the free energy cost of 
creating them - their interface free energy- can be made 
arbitrarily small. 

There exists simulational evidence as to the validity 
of the conclusion that the IRSB state does not exist in 
the 10-state Potts spin glass IQJ. In the infinite-range 
limit when all the Potts spins are coupled by random in- 
teractions, mean-field theory becomes exact and is the 
same as that for the p-spin model. The simulations of 
it were consistent with the existence of both a dynami- 
cal transition Td and a transition to a IRSB state below 
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Tc- However, in three dimensions using a Gaussian dis- 
tribution of the couplings, there was no evidence of a 
dynamical transition Tj^, nor of a transition to a IRSB 
state. 

The complete absence in three dimensions of the dy- 
namical behavior found at mean-field level in the Potts 
glass is perhaps surprising given that in supercooled liq- 
uids there is ample evidence that mode-coupling works 
well. Furthermore we have recently been able to derive 
the replicated form of the p-spin functional starting di- 
rectly from the liquid state However, this derivation 
of thep-spin functional provides the answer to the puzzle: 
Its variables are quantities which are not directly related 
to the densities whose correlations functions are studied 
in mode-coupling theory. The derived replicated p-spin 
functional is useful though when studying whether there 
might be a genuine glass transition. It can be further 
transformed to the replicated functional which describes 
an Ising spin glass in a field ■ Then the issue of 

whether there is or is not a genuine thermodynamic glass 
transition is the same as whether there is a transition of 
an Ising spin glass in a magnetic field. But whether there 
is or is not a transition to a glass state, that state cannot 
be the IRSB state as in the mean-field p-spin glass model, 
since such a state does not exist in finite dimensions. 

Numerical studies of p-spin models outside the mean- 
field limit use Hamiltonians which are of the following 
form, which is appropriate for the case p = 3 [l^ . Each 
site is occupied by two Ising spins, Ui and r^, and the 
Hamiltonian is 

+ jJV,a,r, + j(;V,a,T,) , (2) 

where < ij > are nearest-neighbor pairs, and the cou- 
plings Jj'-"'' are chosen independently from a Gaussian 
distribution with zero mean and width J. Note that 
when the signs of all the spins are reversed, the sign of 
the Hamiltonian changes, indicating the violation of time 
reversal symmetry. For models which have p even, such 
as that for p = 4 in Ref . ^3 1 the Hamiltonian does not 
change sign when all the spins are flipped. 

Starting from a p-spin Hamiltonian such as that in 
Eq. (O it is possible to obtain the 'generic' form of the 
replicated p-spin functional 0, 0, 0| by averaging 
over the couplings: 
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where qaf3 is a symmetric matrix with Qaa — 0. The 



coeflicients t, wi, W2 and y are arbitrary positive param- 
eters. This generic functional has the same relationship 
to the finite dimensional p-spin glass models as (say) the 
Landau-Ginzburg functional has to the Ising model. It 
represents a p-spin model with odd values of p as the 
cubic term with coefficient W2 corresponds to a term 
which breaks time-reversal invariance. (Terms in which 
a replica index appears an odd number of times always 
break time-reversal invariance) . The odd p models would 
seem to be the relevant ones for applications to struc- 
tural glasses which have no identifiable symmetries. The 
replica indices go from 1 to n, and Qaa — 0. 

In order to calculate interface free energies between 
systems which differ as to whether they have periodic 
or anti-periodic boundary conditions, we will follow the 
procedure used in Ref. 9] . This requires one to replicate 
the system with periodic boundary conditions n times 
and the system with anti-periodic boundary conditions 
m times, where as usual we will take n and m to zero at 
the end of the calculation. 

Expanding the replicated partition function in powers 
of m and n, we have 
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{n + to) i5F 



(n + to)^ 



P.AP 



(4) 



where AF^ = Fp — Fp = F\p — Fap is the (mean 
square) sample-to-sample fluctuation of the free energy, 
the same for both sets of boundary conditions P or AP, 
and F = = T^p- Thus AFp^ap = 0. To find the 
variance of the interface free energy, AFp , we expand 
out hiZpZ™p to second order in the numbers of repli- 
cas, n and to, separate out the pieces involving the total 
number of replicas n -\- m, and take the remaining piece, 
which is proportional to nm. 

The replica indices in the generalization of Eq. JSJ go 
a, /?, 7 = 1, 2, • • • , n, n -I- 1, • • • , n + m. The order pa- 
rameter q divides naturally into blocks of size n and to. 
From now on, Greek indices will label the first block, Ro- 
man ones the second block, so, for example, qai, means 
a € [^,n\ and i £ [n+l,n + m], and refers to the respec- 
tive entry in the off-diagonal, or mixed, sector. 

The boundary condition which results from the flipping 
the sign of of the bonds in the plane perpendicular to the 
z-axis is that q must be periodic in the Greek and Roman 
sectors, and anti-periodic in the mixed sectors: 



qapiz) = qapiz + L) 

qtj (z) = qij {z + L) 

qai{z) = -qai{z + L). 



(5) 



In the directions perpendicular to the z-axis, which we 
will take to be all of length M, the system is periodic. 
At mean-field level, the stable solution for InZpZ'^p is 
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given by 



lnZ|iZJp = /3H,.cp{g^P}, (6) 



where q^^ denotes the solution of the Euler-Lagrange 
equations: 



dz^ 



QaiQifS 

7=1 

W2qlp + yqlp/S, 



i=l 



-- tqij -wi^ qikqk] -wi^ Qi-yq-yj 

k=l 7=1 

W2ql +yqfJi, 



(7) 



= tqai -wi^ qaiq-ris - '^i^ qakqki 

7=1 fc=l 

W2qli + yqlt/i- 



Notice that there is a solution of these equations which 
satisfies the boundary conditions in which the mixed 
overlaps qai = and where in the all Greek or all Ro- 
man sectors qap and qtj, the solutions are z independent. 
We believe that this is the physically relevant solution as 
it is stable. Thus at mean-field level the free energy of 
the n + m replicated system with the boundary condi- 
tions breaks up into the sum of the free energies of the n 
and m replicated systems without boundary conditions, 
so that there can be no term of order nm in Eq. Q , so 
that the interface energy vanishes to this order. 

In the high-temperature phase the solution to the sta- 
tionarity equations is qap = qij — qai = 0. A spa- 
tially constant IRSB solution appears discontinuously at 
a temperature tK — 2w'^/{3y) where to simplify the al- 
gebra the following notational changes have been made: 
wi = 1, 'W2 = ^ + w with u; > 0. (For the specific model 
of Eq. Q ui is actually negative 15]). There is, for 
temperatures t < Ik, a IRSB solution in which qap is 
zero everywhere except in n/m boxes of size m along the 
leading diagonal, where it takes the value qi. qij simi- 
larly is zero everywhere except in m/rh boxes of size rh. 
qai = 0. For a review of IRSB see Ref. Right at Ik, 
on letting n,m go to zero, qi = 2w/y and fh — 1. For 
temperatures t < t^, qi increases and m decreases. 

To get the leading non-vanishing contribution to the 
interface free energy we need to go to one-loop order. The 
first correction is due to Gaussian fluctuations around the 
stationary IRSB solution: 

= /3Kep{g^P}+i^/(fc2), (8) 



In ZpZ^p 



fc is a d-dimensional wave vector, are the eigenval- 
ues of the Hessian, evaluated at the stationary point, 
and their degeneracies are d^. The eigenvalues and de- 
generacies are the same as for a system of size n + m 
without boundary condition changes (because the sta- 
tionary point is the same), only the nature of the k- 
vectors changes for the terms involving eigenvalues whose 
corresponding eigenvectors / are nonzero exclusively in 
the mixed sector (i.e. fa/s — fij — 0): the wave vectors 
have to respect the imposed boundary conditions, which 
implies k = {2nin/M, 2nd-in/M, {2nd + 1)71"/^) 
(with Hi € Z) in the mixed sector as opposed to A; = 
(2rti7r/M, . . . , 2nd-iTT/M, 2?i<j7r/L) in the Greek or Ro- 
man sectors. 

While it is possible to compute the eigenvalues and 
their degeneracies d,,, it is much easier to follow the pro- 
cedures of Ref. |e113- One computes dl/d{k'^) rather 
than / itself. This helps because it is possible to express 
the terms in this derivative of order , m? and nm di- 
rectly in terms of a single propagator of the mixed sec- 
tor calculated for the system without boundary condition 
changes: 



dlAP 

a(fc2) 



~ nmGn 



(10) 



where the prefactor nm reflects the number of eigenvec- 
tors in the mixed sector. The subscript AP indicates 
that its argument k must have values appropriate to the 
mixed sector. Defining the integral 



(11) 



J(fc) :== / d(fc^)a 



and obtaining the constant of integration via Ref. 
enables one to write the terms of order n^, m^ and nm 
in / as 



+ m^ ^ ,,2 



Jp{k )+nmJA_p{k ) 



(n + m)' 



Jp(k^) + nm{JAp{k^) - Jpik^))- 



The subscripts P and AP on J mean that J must be 
taken as when the argument is not of the required type, 
i.e. periodic or anti-periodic. 
Comparison with Eq. Q shows 
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where 



/(fc2) = ^d^ln(fc2 + A, 



where the subscripts on the sums indicate the nature of 
the allowed A:- vectors, as made explicit in the second part 
(9) of the equation where the z component of the fc-vector 
has been split off, leaving the d — 1-dimensional wave 
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vector I. The sum over the z component can be extended 
to ±00, introducing only exponentially small errors for 
large L. 

The calculation of Gn^ni is readily accomplished using 
the methods of Ref. [l9j. Right at the transition i^, 
the propagator is very simple because to = 1: Gai,ai = 
+ t), so J(fc2) = ln(fc2 + t). We shall for simplicity 
work in the limit where the sums over the first d — 1 
directions can be converted to integrals (renaming the 
integration variable q). Then 



d-l / " 1 



P-^FI,, - / (2.) 
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+ t)-ln(g^ + M!^+i)). 



The sum over r can be done analytically when the right- 
hand side becomes 

ij^^^M"-^ dqq''-^\ncoth[^/q^tL/2]. 

Sd is the surface area of a d-dimensional unit sphere. In 
the limit when L ^ ^, where — is the correlation 

length, the integral is simple to evaluate: 



P.AP 
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This result is valid in the limits M ^ L ^ By going 
back to Eq. (|12|l . however, any desired aspect ratio can 
be studied for any value of L/£^. 

While we have produced an explicit expression for the 
interface energy SF at one particular temperature, tx, at 
which the system jumps into the IRSB state, the same 
methods can be used at any temperature. The origin 
of the exponentially small free energy is the presence of 
positive non-zero eigenvalues in the Hessian matrix. The 
stationary point needs to be marginally stable in order 
to get interface free energies of the form SF ~ L^, with 
^ > 0, which is what happens for the Ising spin glass 
0. In both the temperature regions above and below 
Ik, the eigenvalues are positive so the interface free en- 
ergy in p-spin models is at all temperatures exponentially 
small. This conclusion will also be true of the IFRSB 
state which exists for some generalizations of the p-spin 
model, where the order parameter q{x) has both a jump, 
characteristic of IRSB, but in addition varies with x, like 
in full replica symmetry breaking. It has also only non- 
zero Hessian eigenvalues [isj . 

The methods used here can also be used to calculate 
the interface free energy appropriate to the dynamical 



questions studied in Ref. 0; 13 • To locate the dynamical 
transition one imposes the condition that the stationary 
point solution q^^ has a Hessian with a null eigenvalue 
|l8l | . This will guarantee that the associated interface free 
energy will not be exponentially small. The calculation 
is very similar to that of Ref. 



liy sma. 
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